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ABSTRACT 



Context. Owing to its unusual location and its isolation, the nature of the high Galactic latitude 09.5 Vp object HD 93521 is still 
uncertain. 

Aims. We have collected X-ray and optical observations to characterize the star and its surroundings. 

Methods. X-ray images and spectra are analyzed to search for traces of a recent star formation event around HD 93521 and to search 
for the signature of a possible compact companion. Optical echelle spectra are analysed with plane-parallel model atmosphere codes, 
assuming either a spherical star or a gravity darkened rotationally flattened star, to infer the effective temperature and surface gravity, 
and to derive the He, C, N and O abundances of HD 93521. 

Results. The X-ray images reveal no traces of a population of young low-mass stars coeval with HD 93521. The X-ray spectrum of 
HD 93521 is consistent with a normal late 0-type star although with subsolar metallicity. No trace of a compact companion is found 
in the X-ray data. In the optical spectrum. He and N are found to be overabundant, in line with the effect of rotational mixing in this 
very fast rotator, whilst C and O are subsolar. A critical comparison with the properties of subdwarf OB stars, indicates that, despite 
some apparent similarities, HD 93521 does not belong to this category. 

Conclusions. Despite some ambiguities on the runaway status of the star, the most likely explanation is that HD 93521 is a Population 1 
massive O-type star that was ejected from the Galactic plane either through dynamical interactions or a result of a supernova event in 
a binary system. 

Key words. Stars: early-type - Stars: fundamental parameters - Stars: massive - Stars: individual: HD 93521 - X-rays: stars 



1. Introduction 

O-type stars, the hottest and most massive main-sequence stars 
of Population I, are seldom found in isolation. They are prefer- 
entially found in open clusters and OB associations, which are 
therefore most likely to be their birth places (e.g. de Wit et al. 
12005 ). The few O-type stars that are not directly located inside 
a cluster or association are usually runaways and are hence be- 
lieved to have been ejected from their parent cluster either by the 
kick of a supernova explosion in a binary system or as a result 
of dynamical interactions in a dense cluster core. 

The high Galactic latitude O-type star HD 93521 {hi = 
183.14°, bii - 62.15°), located very far away from any known 
site of recent star formation, challenges this general picture. 
While the optical spectrum of HD 93521 leads to an 09.5 Vp 
classification, the nature of this star has been subject to debate 
over many years (see e.g. Ebbets & Savage [T982I Irvine [T9891 
Howarth & Smith 2001) and the case is still not settled. 

HD 93521 has one of the largest rotational velocities known 
among Galactic O-stars (390km s ', Rauw et al. GOOSb . The stel- 
lar wind has an apparently low terminal velocity and is likely 
heavily distorted into a Be-like decretion disk wind (Howarth 
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& Reid (19931, Bjorkman et al. [19941 Massa [T995] l. In the opti- 
cal domain, the wind produces emission features in the wings of 
the Ho- line, although they are far less prominent than in gen- 
uine Oe stars. In addition, HD 93521 displays bi-periodic (1.75 
and 2.89 hr) absorption line profile variability that is commonly 
interpreted as the signature of two non-radial pulsation (NRP) 
modes with / a; 8 + 1 and Z w 4 + l(see Rauw et al. 120081 
and references therein). An alternative explanation for the op- 
tical line profile variations, that cannot totally be ruled out a pri- 
ori, would be the effect of a compact companion (formed in the 
supernova explosion that ejected the system from the Galactic 
plane) orbiting the O-star and accreting material from its rota- 
tionally flattened wind. Indeed, in the case of Be/X-ray binaries, 
the decretion disk can be truncated by resonance with the or- 
bit of the neutron star companion (Okazaki & Negueruela l2001l) 
and higher order resonance could possibly trigger periodic den- 
sity waves in the disk that would contaminate the photospheric 
absoption lines via variable residual emission. 



To help clarify the nature of HD 93521, we have collected 
new optical and, for the first time. X-ray observations. Our obser- 
vations are presented in Sect.|2] Section[3]is devoted to the study 
of the surroundings of HD 93521 as seen in our X-ray images. 
The X-ray spectrum of HD 93521 is analysed in Sect.lH whilst 
the optical spectra are studied in Sect.|5] using plane-parallel at- 
mosphere models assuming either spherical or rotationally flat- 
tened geometries for the star. The results are discussed in Sect.|6] 
and Sect.|7]presents our conclusions. 
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2. Observations 

2.1. X-ray data 

A 40ksec X-ray observation was obtained on 2 - 3 November 
2009 with the XMM-Newton satellite (Jansen et al. 120011) . To 
reject optical and UV photons, the EPIC cameras (Turner et al. 
120011 Striider et al. 120011) were used with the thick filter. 

The raw data were processed with the SAS software ver- 
sion 10.0. The end of the observation was affected by the raise 
of a background flare. We discarded this part of the observa- 
tion, ending up with an effective exposure time of 38 ksec for 
EPIC-MOSl, 37 ksec for EPIC-MOS2 and 30 ksec for EPIC- 
pn. Images were built over soft (0.5 - 1.0 keV), medium (1.0 
- 2.0keV) and hard (2.0 - S.OkeV) energy bands with a pixel 
size of 2.5"x2.5". These images were exposure corrected and 
combined into an energy-coded three-colour image of the field 
of view (Fig.[l]). 




Fig. 1. Energy-coded three-colour image built from our XMM- 
Newton observation of HD 93521 (the reddish source near the 
center of the image). The red, green and blue colours correspond 
to the soft, medium and hard energy bands used throughout this 
paper (see text). The individual EPIC images were exposure cor- 
rected before they were combined. A colour version of the figure 
can be found in the electronic version of the journal. 



2.2. Optical spectroscopy 

A series of twenty four optical spectra were obtained over 
four nights in April - May 2004 with the ELODIE echelle 
spectrograph (Baranne et al. 1 19961 1 at the 1 .93 m telescope at 
Observatoire de Haute Provence (OHP). ELODIE was an echelle 
spectrograph with a resolving power ofR^ 42000 over the 3850 
- 6800 A wavelength domain. Our initial intention was to use 
these data for our study of the line profile variability (Rauw et 
al. 120081 ). However, due to the very poor weather conditions dur- 



ing this run, individual spectra lacked a sufficient signal-to-noise 
ratio and we thus decided to combine all the spectra into a sin- 
gle high-quality spectrum. This average spectrum, with a signal 
to noise ratio (S/N) of 390 at 5500 A, was obtained from data 
collected at many different phases of the line profile variation 
cycles and should hence be relatively free of any signatures of 
the latter. 

Another echelle spectrum was obtained in November 2010 
with the SOPHIE spectrograph (Perruchot et al. 120081 ) also at 
the 1.93 m telescope at OHP. SOPHIE covers the wavelength 
domain from 3870 to 6940 A. The spectrum was reduced with 
IRAF. This spectrum is of higher resolution than the ELODIE 
ones with a resolving power, R ^ 74 000, as measured from cal- 
ibration lamps. The S/N of this spectrum is ^ 310 at 5500 A. 

2.3. Ha imaging 

Since there are apparently no Ho- images available for the re- 
gion around HD 93521, we collected some images with the help 
of two amateur astronomers, Gaston Dessy and Joel Bavais. 
For this purpose, G. Dessy used a TMB-92 9.2 cm refractor 
with a focal length of 40.8 cm equipped with an Atik 16IC 
CCD (659 X 494 pixels), and for a second run with an Atik 
4000M CCD (2048 x 2048 pixels). Both CCDs have pixel sizes 
of 7.4 X 7.4jum^, corresponding to 3.7"x3.7"on the sky. J. 
Bavais used a SkyWatcher 80/600 ED refractor equipped with 
an Atik 31 4L CCD (1392 x 1040 pixels of 6.45 x 6.45 i^im^ corre- 
sponding to 2.2"x2.2"on the sky). In all cases, commercial Ha 
(70 A bandwidth) and V-band filters manufactured by Baader 
Planetarium were used. For each set-up, integration times were 
1 hour in Ha (split into 12 exposures of 5 min each) and 15 min 
in V (split into 5 exposures of 3 min each, except for the last 
run where no V-band images were taken). All data were pro- 
cessed with the MIDAS software developed at ESO. The surface 
brightness of the sky measured on our images amounts to about 
17.8 mag arcsec"^ in V. After subtracting a flat sky background, 
the residual 1-cr fluctuations in the V band images correspond to 
a surface brightness around 22.0 mag arcsec"^. The Ha images 
are about 1 .5 mag shallower than the V-band data. No trace of 
a nebular emission was found on the images before or after sky 
subtraction. 

3. The surroundings of IHID 93521 

The position of HD 93521 far away from the Galactic plane and 
from any open cluster naturally raises the question of the origin 
of this star. There are currently only a few known examples of 
massive stars that have likely formed in isolation. For instance, 
de Wit et al. (2005 ) conclude that 4 + 2% of the Galactic 0-stars 
could have formed outside cluster^]. 

Models for massive star formation based on the competitive 
accretion scenario imply that a high-mass star must be associated 
with a population of low-mass stars (Bonnell et al. 20040. Core 
accretion models of massive stars on the contrary allow rather 

' In the same context, Bressert et al. ( 120 12t recently identified a sam- 
ple of massive stars, earlier than spectral type 07, within a projected 
distance of 125 pc around the 30 Dor complex, that are not associated 
with any cluster and are thus candidates for massive stars formed in iso- 
lation. It has to be stressed though that the case of HD 93521 requires 
a far more extreme form of isolation than considered by Bressert et al. 
l f20T2b . 

- Although dynamical interaction could lead the sparsest clusters to 
disperse rather quickly, on time scales of a few Myr. 
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isolated O-stars to form (Krumholz et al. l2009l l. Isolated O-type 
stars that are still in their formation region are surrounded by 
residual gas that can be seen through Ha observations (Lamb 
et al. 120101 Selier et al. 1201 11 1 and some of them are found in 
sparse clusters with less than ten lower-mass companions (Lamb 
et al. 120101) . Monte Carlo simulations by Lamb et al. (1201 Oi l in- 
dicate that the existence of such sparse clusters is more in favour 
of the core accretion models and suggest that 'clusters are built 
stochastically by randomly sampling stars from a universal ini- 
tial mass function (IMF)'. 



0.01 



0.001 



0.000 




0.5 1 

Energy (kcV) 

Fig. 2. Best-fit absorbed power-law model (see Table[T]) of the 
X-ray spectrum of the blazar BZQ J1049-I-3737. The black and 
blue data points refer to the pn and MOS2 data, respectively. The 
lower panel shows the contribution of the various energy bins to 
the;i'^ of the fit times the sign of the difference between the data 
and the model. 



X-ray observations can help us uncover low-mass, optically 
dim members of very young open clusters. Indeed, X-ray im- 
ages of massive stars in dense open clusters usually reveal a 
rather rich population of X-ray bright, low-mass pre-main se- 
quence stars (e.g. Damiani et al. 2006, Sana et al. 2007). This 
statement also holds to some extent for rather scarce open clus- 
ters such as NGC 6383 (Rauw et al. '2003^ and to less densely 
populated OB associations (e.g. Rauw 2011). The brightness of 
very young low-mass stars in the X-ray domain makes this the 
ideal energy range to distinguish between foreground or back- 
ground field stars and genuine cluster members coeval with the 
massive stars. This technique was successfully applied to reveal 
the population of low-mass pre-main sequence (PMS) stars asso- 
ciated with two rather isolated BO stars in the diffuse H n regions 
S 255 and S 257 (Mucciarelli et a HIOTlT l. 

By analogy with the two BO stars investigated by Mucciarelli 
et al. (201 1), one expects about 300 low-mass pre-main sequence 
stars associated with HD 93521 if this stellar population follows 
a standard IMlfl We can now estimate what fraction of these 
low-mass objects would be detected with our X-ray observation 



and compare this to the actual number of sources in the field of 
view. 

The three-colour image of the field of view around HD 93521 
(Fig-Hi' does not reveal any obvious concentration of X-ray 
sources around HD 93521. Applying the SAS source detection 
routines with a significance thresholcQ of 10 yields a total of 
55 sources. After visual inspection, 10 sources were removed 
from the list because their detection was affected by gaps be- 
tween CCDs or because they were at the edge of the field of 
view. Two more faint sources that were apparently missed by 
the detection algorithm were added to our final list of 47 X-ray 
sources. Assuming an optically thin thermal plasma model with 
a typical temperature of 1 ke V, our detection limit corresponds to 
an observed flux of about 3 x 10"'^ ergcm"^ s"'. At an adopted 
distance of HD 93521 of 1.2 kpc (see below), this flux would cor- 
respond to an X-ray luminosity of 5 x 10^^ erg s ' . If we compare 
this number with the cumulative X-ray luminosity functions of 
the Orion Nebular Cluster (Preibisch & Feigelson ,2005i) , we es- 
timate that we should detect about 90, 85 and 50% of the PMS 
stars in the mass ranges 0.9 -1.2, 0.5 - 0.9 and 0.1- 0.5 M© re- 
spectively. Therefore, we would expect to observe of order 100 

- 200 X-ray sources in the vicinity of HD 93521. This is clearly 
not the case here. 

To assess the impact of background X-ray sources, we have 
repeated the detection chain on two different energy bands 0.5 

- 2.0keV and 2.0 - 8.0keV. In the 0.5 - 2.0keV band, 42 
sources are detected with a limiting flux of order 3 - 4 x 
10^'^ergcm"^ s"'. From the logA^ - log 5 relation of Giacconi 
et al. (1200 11 1 derived from the Chandra Deep Field South, we 
would expect about 50 detections of extragalactic point sources 
in this energy band and with this flux limit. This result implies 
that the bulk of the detected sources must be extragalactic back- 
ground objects, unrelated to HD 9352 1 . Very similar conclusions 
are reached with the results for the 2.0 - 8.0 keV band. This 
conclusion is further supported by the fact that only five X-ray 
sources in our field of view have a 2MASS near-IR counterpart 
within a correlation radius of 4 arcsec. Out of these five sources, 
four are among the brightest X-ray sources: HD 93521 itself, the 
blazar [MGL2009]BZQ J1049-H3737 (Massaro et al. 2009), the 
stai- BDh-38°2183 (with a SIMBAD spectral type F8), and an 
unknown source at 10:48:14.9 -1-37:23:24 (hereafter source 4). 

We have analysed the EPIC spectra of the three brightest 
sources in addition to HD 93521. The results are listed in Table[T] 
The blazar and the unknown source are best fitted with a power- 
law model (see Fig.|2]i, suggesting that the latter is most proba- 
bly also related to an AGN. The spectrum of the late-type star is 
well represented by a single temperature optically thin thermal 
plasma model. 

We thus conclude that our XMM-Newton data do not reveal 
any evidence for a lower mass stellar population that could be 
associated with the formation of the O-type star HD 93521. This 
is in line with the lack of nebular emission from residual gas in 
our Ha images (see Sect. l2.3l l. 



4. The EPIC spectrum of HD 93521 

HD 93521 is detected as a moderate X-ray source with net 
count rates of 7.3 x 10"^ 8.2 x 10"^ and 3.6 x 10"^ ctss"' 



^ This estimate is in line witii the relation between the total cluster 
mass and the mass of its most massive member as proposed by Weidner 
et al. ( 1201 Oi l. Indeed, applying this relation to HD 93521, one would 
expect a cluster with a total mass in the range 100 - 150 M©. 



This implies a probability of < e"'" that a random Poissonian 
fluctuation could have caused the observed source counts inside the 
5x5 pixel^ detection cell. Over the entire EPIC field of view, the ex- 
pected number of spurious detections amounts to about 2% of the total 
number of detected sources. 
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Table 1. Spectral properties of the brightest X-ray sources other than HD 93521. 



Source 


Model 


(10^2 cm-2) 


kT 
(keV) 


r 


fx {0.5- lOkeV) 
(erg cvcT^ s"') 


/^""^ (0.5 - lOkeV) 

(ergcm^^ s"') 


xi 


d.o.f. 


BD+38°2183 


wabs*apec 


< 0.020 






1.11 X 10^'^ 


1.11 X 10^'' 


0.90 


11 


BZQJ1049+3737 


wabs*power 


024+"-" 




1.67M" 


2.12 X 10-'^ 


2.18 X 10^''' 


1.05 


65 


source 4 


wabs'-power 


0.045+;;j[^ 




2.02+;|3 


2.27 X 


2.43 X 10-'3 


1.54 


73 



for EPIC-MOSl, MOS2 and pn, respectively. The EPIC spec- 
tra of HD 93521 were analysed using the xspec software (ver- 
sion 12.6.0, Arnaud [l996l l. Unless stated otherwise, the inter- 
stellar neutral hydrogen (Hi + H2) column density was set to 
1 .3 X 10^° cm"2 (Bohlin et al. |1978l l. To account for the possible 
presence of an additional circumstellar absorption column due to 
the partially ionized stellar wind, we further included an ionized 
wind absorption model (Naze et al. I2004I I with the wind column 
density treated as a free parameter. 

We have tested a variety of models that are potentially ade- 
quate depending on the actual nature of HD 93521. 

Medium resolution X-ray spectra of normal O-type stars are 
usually well fitted with absorbed optically thin thermal plasma 
models (see e.g. Naze 2009). A model with a single plasma com- 
ponent does not fit the EPIC spectra of HD 93521 simultane- 
ously at low and high energies. We have therefore tested two- 
temperature apec models (Smith et al. I200II I. either with solar 
metallicity or with the global metallicity of the plasma taken as 
a free parameter (see models [1] and [2] in Table|2|i. In the latter 
case, the best fit is achieved with a metallicity near 0.22. This 
fit provides a significant improvement over the solar metallicity 
solution. 

Using the He and CNO abundances derived from our best fit 
of the optical spectrum (see Sect.|5]l, we have then tested a two 
temperature apec model with variable abundances (vapec). The 
He, C, N and O abundances were fixed at respectively 2. 14, 0. 19, 
1.51 and 0.30 times solar. As a first step, all other elements were 
taken to have solar abundances. This led to a relatively poor fit 
with^J = 1.52. Apart from C, N and O, Fe and Ne have strong 
lines in the X-ray domain covered by our EPIC data. Therefore, 
we allowed the abundances of these two elements to vary in the 
fit. This results in an improved fit with Ne and Fe abundances of 
respectively 0.28;^ ?^and 0.43+jg'* relative to solar (see Fig.[3]and 
model [3] in TablelJ]). Assuming that Ne and Fe have the same 
abundances as O (i.e. 0.30 solar), which is our best indicator of 
the metallicity as the O abundance is only little affected by the 
CNO process, yields an equivalent fit (see model [4] in Table|2]). 

Whatever the metallicity, we find that the spectrum can be 
represented by two plasma components of temperatures near 
0.3 and 3.0keV, though the hotter component is poorly con- 
strained (especially towards higher temperatures). We further 
find that a significant wind absorption is needed to provide 
a reasonable fit and that the interstellar absorption corrected 
flux corresponding to thermal plasma models is of order 6.0 - 
6.5x10"'"* erg cm"^ s"' . Assuming a distance of 1.2kpc, this flux 
corresponds to an X-ray luminosity of 1.0 - 1.1 xlO^' ergs"', 
which translates into an Lx/L^oi ratio of 8.7 - 9.4 xlO"**. This 
number is towards the lower end of, but compatible with, the 
range of Lx/^^boi values measured for normal O-type stars with 
XMM-Newton (Naze l2009b . 

A plasma component as hot as 3 keV is somewhat surpris- 
ing for intrinsic X-ray emission by a single ordinary O-type star 
This is especially true if the star has indeed a low equatorial wind 
velocity, of order Vk, - 400 kms"', as derived by Howarth et al. 



d 1 9971 although the polar wind component is likely much faster, 
Voo = 2000 km s ", see Howarth & Reid T993). Naze ( [20091 ) 
found that many O-type stars feature a faint hot plasma con- 
tribution around 2keV, but only a few objects have a plasma 
that is actually as hot as in the case of HD 93521. Over recent 
years, a new category of intrinsic hard X-ray sources associated 
with early-type stars was identified, the so-called j Cas analogs 
(e.g. Smith et al. 12004 Rakowski et al. 120061 1. These are late 
O or early B-type emission-line stars which have X-ray spec- 
tra dominated by plasma components with kT up to 12keV, lo- 
cated rather close to the Be star and its disk. However, unlike 
HD 93521, these objects are highly overluminous (by a factor 
^ 40, Rakowski et al. ,2006) with respect to the typical Lx/L\,o\ 
ratio of OB stars. Therefore, HD 93521 has rather normal X-ray 
properties for an O-type star, except for the somewhat higher 
secondary plasma temperature. 

As an alternative, we have tested the possibility that the hard 
emission might be non-thermal. Such a non-thermal emission 
could arise either from an accretion process in a binary sys- 
tem featuring a compact companion or from inverse Compton 
scattering of photospheric UV photons, if the wind of HD 93521 
contains a population of relativistic electrons. The fit is of sim- 
ilar quality to the best fit purely thermal model (see model [5] 
in Table|2|i. The softer emission still requires a thermal plasma 
component at about 0.3 keV. For this kind of model, we have also 
attempted to let the metallicity of the apec component vary dur- 
ing the fit. There is no improvement of the fit quality and, apart 
from the metallicity which goes to 0.29;^ Z© and the normal- 
ization parameter, there are no changes in the best fit parameters. 
Finally, we also tested a model with the individual abundances 
of He, C, N and O set to the values derived from our best fit of 
the optical spectrum and letting the Ne and Fe abundances vary. 
Again, there is no improvement of the fit and the abundances of 
Ne and Fe converge to values very similar to those obtained with 
model [3]. 

Finally, we have tested a model consisting of a black-body 
and a power law. This kind of model was found to adequatly de- 
scribe the EPIC spectra of the sdO + compact companion system 
HD 49798 (Tiengo et al. 2004) . The results are Hsted as model 
[6] in Table|2| If we keep the interstellar column fixed, the best 
fit is of lower quality than what can be achieved with thermal 
plasma models. We stress that the best-fit black-body component 
of HD 93521 would be much hotter (kT = 0.15 vs. 0.034 keV) 
and the power law significantly steeper (F = 1.4 vs. 2.0) than in 
the case of HD 49798 (Tiengo et al. 120041 ). A significantly better 
fit (xl = 1 .08) could be obtained with the black-body + power 
law model if the interstellar column density would be treated as 
a free parameter. However, in this case, the fitted ISM column 
density would reach 0.14 x 10^^ cm"^, i.e. it would exceed the 
observationally determined column density by a factor 10. Such 
a large discrepancy seems rather unlikely. 

To further test the possibility that HD 93521 could be a bi- 
nary system hosting a compact companion (either a white dwarf 
or a neutron star), we have searched the X-ray data for period- 
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Table 2. Spectral fits of HD 9352 1 . 



Model 


log N„ind 

(cm-2) 


kT, 

(keV) 


normi 




kT2 
(keV) 


r2 


norm2 




z 

Zo 


fx 

(10-'* 


fx 

ergcm-^ s-') 


Xi (d.o.f.) 


111 

J 


21.63! 


08 
09 


27+ 


rO 89+" '''^ 10" 


4 


3 01^^-^^ 

-'•"'^-0.96 




CO 18+"' ^ 10- 


-4 


1 (fixed) 


6.20 


6.52 


1.49 (65) 


[2] 


21.54! 


12 
21 


0.28!;°2 


(2.65!i;^) 10" 


4 






(o.22!;j;;^) lo- 


-4 


22+ 

"•^ -.09 


5.89 


6.19 


1.34 (64) 


[3] 


21.37! 


18 
23 




(1.15!!J;g)10- 


4 


2 96+'° ''^ 

^■^"-1.35 




(o.i4!„"^,) lo- 


-4 


model 


5.57 


6.02 


1.32 (63) 


[4] 


21.26! 


17 
33 


o.30!;;;i 


(1.08!jJ;^^) 10- 


4 


4 03+" ' 




co. i2!;j;;!5) lo- 


-4 


model 


5.72 


6.00 


1.31 (65) 


[5] 


21.54! 


09 
11 


0.28!;!| 


(0.62!°;?3) 10- 


4 




2 54+-2'' 


co. io!:Sj^) 10- 


4 


1 (fixed) 


5.97 


6.27 


1.35 (65) 


[6] 






152+""** 


(8.5!!J;?) 10-" 






1 44+36 


(3.0![;^) 10-' 






6.15 


6.41 


1.46 (66) 



Notes. All fits were performed with the interstellar neutral hydrogen column density set to 1.3 x 10^" cm-^. Models [1] and [2] correspond to 
wabs*wind*apec (2T) with solar and free metallicity respectively. Models [3] and [4] correspond to wabs -'wind*vapec (2T) with the abundances 
of He, C, N and O set to 2.14, 0.19, 1.51 and 0.30 times solar. In model [3], the abundances of Ne and Fe are free (see text), whilst they are fixed 
at 0.30 solar in model [4]. Models [5] and [6] correspond to wabs"'wind*(apec+power) and wabs"'(bbody+power) respectively. The fluxes are 
evaluated over the energy range 0.5 - lOkeV. The normalization of apec models are defined as j dV where d is the distance in cm, 
«(, and iiH are the electron and proton densities and V is the volume of the emitting plasma. The normalization parameters of the power law and 
black-body component are respectively the number of photons keV-' cm-^ s-' at 1 keV and the source luminosity in units 10^' erg s-' divided by 
the square of the distance in units of lOkpc. 
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Fig. 3. Same as Fig.|2] but for the best-fit model [3] (see Table|2]i 
of the X-ray spectrum of HD 93521. The black, green and blue 
data points refer to the pn, MOS 1 and MOS2 data respectively. 



icities that could be related to X-ray pulses from the compact 
companion accreting material from the 09.5 component. For 
this purpose, we have extracted event lists from the source re- 
gion of HD 93521 (containing photons from the source and the 
backgrouncQ) for each of the three EPIC instruments. The ar- 
rival times of these photons were corrected into the barycentric 
frame of reference. We then folded the arrival times of all the 
photons with a trial period and built a histogram of the resulting 
arrival phases (adopting phase bins of 0.05). We finally compute 
the maximum difference between the number of counts in the 
phase bins of the histogram. The whole process was repeated for 
10000 trial periods between 1.0 and lOOOfl 

and the maximum 

difference in the number of counts was plotted as a function of 



' The source accounts for roughly 70% of the photons, the remainder 
are due to the background. 

^ In the full frame modes used in our observation, the pn and MOS 
detectors have a time resolution of 0.0734 and 2.6 s respectively. 



the trial period. If there was a periodic signal in the X-ray data 
of HD 93521, we would expect a peak in this diagram. No such 
peak is found in our data (see Fig.|4]i and we thus conclude that 
the X-ray emission of HD 93521 does not contain any signature 
of short-period pulsations that could hint at the presence of an 
accreting compact companion. 




Fig. 4. Maximum difference between the number of X-ray 
counts recorded per 0.05 phase bin as a function of the trial 
period. The various panels yield the results for the three EPIC 
detectors. The total number of counts (source + background) 
recorded in the source region were 363, 362 and 1 143 respec- 
tively for the MOSl, MOS2 and pn detectors. 
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5. Analysis of the optical spectrum 

5.1. Previous work 

The optical spectrum of HD 93521 has been investigated previ- 
ously with model atmosphere codes. The most sophisticated ap- 
proach so far, was the work of Howarth & Smith (2001 ). These 
authors used non-LTE, hydrostatic, plane-parallel H/He model 
atmospheres to analyse the spectrum of HD 93521, accounting 
for the variation of the local gravity (and hence temperature) 
as a function of stellar latitude at the surface of a rotationally 
distorted star. Howarth & Smith (2001) assumed an inclination 
of 90° and their best-fit model yielded Tj^^ = 38000 + 1500 K, 
r^jj. = 33500 K, log gP = 3.9 + 0.1 and log = 3.5 (in cgs 
units), a helium abundance by number of y = 0.18 + 0.03 and a 
star rotating at 90% of its critical velocity. Here the '/?' and 'e' 
superscripts refer to the parameters respectively at the poles and 
at the equator of the star. 
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Fig. 5. Comparison for HD 93521 between the SOPHIE (black) 
and best-fitting He synthetic line profiles (solid red). The dashed, 
red lines show the line profiles computed for the final parame- 
ters derived from the analysis of the SOPHIE spectrum (Table 
[3]). The light grey-shaded areas delineate the regions where the 
quality of the fit has been evaluated. 



dances of He, C, N and O. These results are used in Sect. l5.3l as 
starting points for a model where we account for the effects of 
gravity darkening on a rotationally flattened star. 

5.2.1. Atmospheric parameters and helium abundance 

The atmospheric parameters (rgff and \ogg) and helium abun- 
dance by number, y - N(lie)/[N(li)+N(lie)], have been esti- 
mated by finding the best match between a set of observed H 
and He line profiles and a grid of rotationally-broadened, syn- 
thetic profiles computed using the non-LTE line-formation code 
DETAIL/SURFACE coupled with LTE Kurucz models (Butler 
& Giddings 1985; Giddings il981j FI Such a hybrid approach 
(LTE atmospheric models, but a full non-LTE treatment for the 
line formation) has been shown to be adequate for late O and 
early B-type stars on the main sequence (Nieva & Przybilla 
120071 Przybilla et aL lIOTTT l. 

The analysis has been independently carried out on the 
SOPHIE and mean ELODIE spectrum. The results are listed sep- 
arately in the following but, as only relatively small differences 
are found, the values eventually adopted will be the averaged 
ones. A SOPHIE spectrum of the naiTow-lined 09 V star 10 Lac 
was continuum normalised and broadened to v sin / = 390 km 
s~', as found for HD 93521 based on Fourier techniques (Rauw 
et al. 2008). The continuum regions were subsequently defined 
and used to normalise the two spectra of HD 93521. 

A set of unblended helium features (6 He i and 3 He n linesj^ 
was chosen after an inspection of a spectral atlas and the 10 Lac 
spectrum which was taken as reference. Three narrow-lined stars 
were analysed to validate the procedures used to derive the atmo- 
spheric parameters and abundances: lOLac (09 V; vsin/ = 25 
km s-i), T Sco (B0.2 V; v sin / = 8 km s ' ) and HD 57682 (09 IV; 
V sin / - 25 km s '). Such a differential approach is expected to 
minimise the systematic errors (see, e.g., Vrancken et al. ,1997] l. 
These stars have their parameters and abundances derived using 
the same codes from classical techniques whereby Tefj is based 
on the ionisation balance of various elements (namely He, C, N 
and Si), log g on the fit of the wings of the Balmer lines and the 
abundances on curve-of-growth techniques (lOLac: this study; 
tSco: Hubrig et al. l2008l HD 57682: Morel IMB- The H and 
He features used may be contaminated by some weak metallic 
lines. The synthetic spectra were computed taking these lines 
into account and using the abundances derived from the curve- 
of-growth analysis. Some portions of the line profiles encom- 
passing metal lines not modelled by DETAIL/SURFACE were 
masked out. For HD 93521, the same metal abundances and mi- 
croturbulence (^ = 10 km s"') as for 10 Lac were adopted. 
Model atmospheres with a helium abundance twice solar were 
used for HD 93521 in accordance with the high abundance found 
(see below). 

An iterative scheme was used to estimate Teff , log g and y. 
The temperature is taken as the value providing the best fit to 
the He i and He ii features with the same weight given to these 
two ions. The gravity is determined by fitting the wings of the 
Balmer lines and y by fitting the Hei features. This procedure 
was repeated until the gravity used to fit the helium lines was 
identical to the value yielded by the Balmer lines, and proved 
for the narrow-lined stars to provide the best match to the refer- 



5.2. 1-D model atmosplieres 

In this section we first analyse the spectrum assuming a spheri- 
cally symmetric star to derive stellar parameters as well as abun- 



' The spectra were convolved using the R0TIN3 
program in the SYNSPEC suite of routines (see 
'http : / /nova .astro . umd . edu/ Synspec4 3/ synspec . html I. 

The relevant lines are Hei ^^4026, 4388, 4471, 4713, 4922, 5016 
and He 114542, 4686, 5412. 
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Table 3. Atmospheric parameters and metal abundances (on the scale in which loge[H] = 12) for the reference stars and for 
HD 93521. 
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Notes. The reference values based on classical methods are shown in boldface. For the standard stars, the second line yields the results of our 
analysis of the actual spectra of these stars using the same approach as for HD 93521, whilst the third line gives the corresponding results for the 
spectra degraded to the v sin / of HD 93521. The solar [N/C] and [N/O] abundance ratios are -0.60 and -0.86, respectively (Asplund et al. i2009t . 
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Fig. 6. Comparison for HD 93521 between the observed (black) and best-fitting synthetic metal line profiles under the assumption 
of a spherical star (red) and a gravity darkened rotationally flattened star (blue). The light grey-shaded areas delineate the regions 
where the quality of the fit has been evaluated. The top panels show the synthetic profiles computed for Tgif = 31 000 K, logg = 3.7 
dex, ^ = 10 km s ' and the abundances obtained using the ELODIE spectrum and assuming a spherical star. 



ence values obtained using classical methods. As already noted 
for HD 93521 (Howarth & Smith 2001), the fit was much better 
for He than for the other Balmer lines. However, the gravity esti- 
mated from He was indistinguishable from the mean value com- 
puted taking all lines into account. This analysis has also been 
repeated after convolving the observed spectra of the reference 
stars with a rotational broadening function with v sin / = 390 km 
s"' . The results are shown in Table[3] For the three standard stars, 
the mean diff'erences between the reference values and the results 
obtained through the synthesis are: Areff = +155+733 K, Alogg 
= -0.09+0.07 dex and Ay = +0.005+0.016. We finally obtain the 



following parameters for HD 93521 assuming a spherical star: 
Teff = 30900+700 K, logg = 3.67+0.12dex andy = 0.18+0.02 
(statistical errors estimated based on the differences for the three 
standard stars between the parameters found and the reference 
ones, as well as the values derived for HD 93521 using the two 
spectra). These results confirm the helium enrichment found by 
Lennon et al. (1991; y = 0.20 + 0.05) and Howarth & Smith 
(2001; y = 0.18 + 0.03). A comparison for HD 93521 between 
the observed and best-fitting He synthetic line profiles is shown 
in FigH] No significant differences are found between the tem- 
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log t (O) 

Fig. 7. Variation of the fit quality in the region 4060-4082 A for 
different combinations of the C and O abundances (colour coded 
as a function of the reduced in the case of the SOPHIE spec- 
trum of HD 93521. The best fit is found for loge(C) = 7.70 dex 
and loge(O) = 8.20 dex. 



perature and helium abundance yielded by the He i singlet and 
triplet transitions (see Najarro et al. 120061) . 

5.2.2. Metal abundances 

This is the first time that the metal abundances are derived for 
this star. Determining its chemical composition could help es- 
tablishing its birth place, as significantly subsolar abundances 
are expected if it was formed in situ far from the Galactic plane. 

The high rotational velocity of HD 93521 makes a selec- 
tion of useful metallic features excessively difficult. As done 
for the hydrogen and helium lines, a number of spectral do- 
mains were carefully chosen and selected on the basis that: (1) 
at most two species should significantly contribute to the blend; 
(2) all the lines with a significant strength should be theoret- 
ically modelled; (3) all the elemental abundances derived for 
the three reference stars (either using the observed spectrum or 
the one broadened to the v sin ; of HD 93521) should reasonably 
match the reference values derived from the curve-of-growth 
analysis (the absolute difference is on average 0.09 dex and at 
most 0.23 dex). Only three spectral regions fulfilled these crite- 
ria: 4060-4082 A (main contributors Cm and On), 4691-4709 
A (main contributor On) and 4995-5011 A (main contributor 
Nil). Unfortunately, the abundances of the other metals (e.g., 
Mg, Si) could not be reliably investigated. The parameters in 
Table [3]have been used, such that any eiTors in Tgff and log g are 
propagated to the derived abundances. The results are shown in 
Table|3] while the observed and synthetic line profiles are shown 
in Fig|6] Figure Q shows an example of the variation of the fit 
quality in the region 4060-4082 A for different combinations of 
the C and O abundances. Taking into account the differences for 
the three reference stars between the abundances found and the 
reference ones, the different values obtained for oxygen from the 
two different spectral domains (there is a good agreement with 
a difference of at most 0.16 dex), and finally the values derived 
for HD 93521 using the two spectra, the statistical errors on both 
the abundances and abundance ratios can be set at the s;0.15 dex 
level. 

Compared to the values derived for the three reference stars 
based on the spectra broadened with v sin / = 390 km s ', the 



abundances of HD 93521 are on average lower by 0.69 dex for 
C, higher by 0.26 dex for N and lower by 0.27 dex for O. A 
nitrogen enrichment at the surface is expected given the high 
helium abundance and is indeed observed: the [N/C] and [N/O] 
abundance ratios are much higher than in 10 Lac and HD 57682, 
and also higher than in r Sco which is the prototype of the class 
of slowly-rotating, main sequence B stars (unexpectedly) dis- 
playing a nitrogen excess (e.g.. Morel et al. 2008). The quan- 
titative behaviour of the [N/C] and [N/O] abundance ratios is 
within the errors what is expected for CNO-cycled material 
dredged up to the surface (Przybilla et al. 120 101) . Determinations 
of the CNO abundances for such fast rotators are very rare in 
the literature: Villamariz et al. (2002J found [N/C]>-i-0.78 and 
[N/O]>-h0.12 dex for HD 191423 (09 III:n; vsin/ = 450 km 
s-'), while Villamariz & Herrero (20051) found [N/C]~+0.48 and 
[N/O]«i-0.35 dex for ^Oph (09.5 Vnn; v sin / = 400 km s '). 

Table 4. Sensitivity of the results to differences in continuum 
placement and choice of microturbulence. For the results with 
a different continuum location, the abundances have been esti- 
mated using the revised parameters. 
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As the abundances of the three metals investigated are af- 
fected at different extent by nuclear CNO-processing, it is diffi- 
cult to estimate the metallicity of HD 93521. However, oxygen is 
expected to be only slightly depleted at these levels of nitrogen 
enrichment and may be taken as a proxy for the global metal con- 
tent. Assuming a depletion of oxygen of ^ 0.1 dex (e.g., Heger 
& Langer 2000 ), there is thus some evidence for a lack of metals 
in HD 93521 compared to solar at the 0.2 dex level. It remains 
to be seen, however, if these low abundances are not an artifact 
of the data treatment: given the lack of true continuum regions, 
there may be a tendency to systematically place the continuum 
too low and therefore to underestimate the strength of the metal 
lines. To investigate this issue, we have redetermined the abun- 
dances of HD 93521 after uniformly shifting the continuum level 
upwards by 0.5%. This leads to an upward revision of T^ff and 
\ogg, as well as abundances that are typically 0.15 dex higher 
(this translates to ^ 0.30 dex for a shift of 1%; Table |4|. We 
therefore conclude that the uncertainties in the metal abundances 
(especially owing to the difficulties in the continuum placement) 
prevent a clear conclusion to be drawn solely based on chem- 
ical arguments regarding the birthplace of HD 93521 (either in 
or far from the Galactic plane). This is especially true when one 
considers that, although the differences are small, the assump- 
tion of a gravity darkened rotationally flattened star leads to sys- 
tematically higher abundances (see Sect. 15. 31 ). It should also be 
kept in mind that the microturbulent velocity is not constrained 
and has been fixed to ^ = 10 km s"'; adopting a lower value 
would result in generally slightly higher abundances (see Table 
|4]i. Nevertheless, our conclusions regarding the He and N-rich 
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status of HD 93521 appear robust against differences in contin- 
uum placement and choice of microturbulence (Table|4ll. 



5.3. Assuming a gravity darkened rotationaiiy fiattened star 

To refine the metal abundances, we have accounted for gravity 
darkening and rotational flattening using an approach similar to 
the one of Howarth & Smith (1200 11 1. 

The stellar surface is assumed to be an equipotential surface 
of the Roche potential model. The latter can be written in spher- 
ical coordinates centred on the centre of a star of mass ot,. 



Q 



G m, o)^ 



— X r sin 0, 
KO) 2 



(1) 



where r is the distance from the centre of the star, 6 is the co- 
latitude, (fi is the longitude and o) the (constant) angular rotation 
velocity. 

The local gravity is then equal to the gradient of the Roche 
potential. We have also take into account the inner radiation pres- 
sure following the approach of Howarth (I1997I I. 



P \4 1 
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i, o-Th = 0.036 m^ kg-' is the 



where F = ^crC^^)' 

Thomson scattering cross section, and cr the Stefan-Boltzmann 
constant. 

The temperature is computed using the von Zeipel ( |1924| | 
theorem 
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where a is traditionally taken equal to 1 for massive stars. 
Recent interferometric observations of rapidly rotating B stars 
(e.g. Kraus et al. 2012), as well as theoretical work (Espinosa 
Lara & Rieutord 2011) suggest however that this might not be 
totally appropriate. Observations are better fitted by a value of 
a » 0.75, whilst the theoretical results suggest a more complex 
law with a dependence of the best-fit a parameter on the actual 
flattening of the star In the following, we have adopted a - 0.75. 
This value of the gravity darkening exponent reduces the temper- 
ature contrast between the pole and the equator compared to the 
traditional von Zeipel law. 

The stellar surface is discretized into 240 x 60 (longitude x 
colatitude) constantly spaced points. In each surface point, the 
local temperature and gravity are known, thus, we can compute 
a local contribution to the spectrum by linear interpolation be- 
tween four spectra of a grid of synthetic spectra computed with 
the non-LTE line-formation code DETAIL/SURFACE coupled 
with the LTE Kurucz models. The spectra of the grid are spaced 
by lOOOK in T^ff and 0.1 dex in \ogg. The appropriate Doppler 
shift is then applied to the local spectra accounting for the rota- 
tional velocity and these spectra are also multiplied by the area of 
the surface element projected along the line of sight towards the 
observer The last corrective factor applied to the local contribu- 
tion to the spectrum is the limb-darkening. The limb-darkening 
coefficient is based on the tabulation of Al-Naimiy (1978) for a 
linear limb-darkening la\\0. The total spectrum is then computed 



' We have compared the synthetic spectra with those simulated us- 
ing more recent tabulations of the limb-darkening coefficients of Claret 
& Bloemen ( 1201 It . The differences between the normalized synthetic 
spectra are of the order 3 x 10"* continuum units, i.e. negligible in view 
of other uncertainties. 



by summing the incremental contributions of each surface point. 
We assume that there is no cross-talk between the different sur- 
face elements. 

This model was then used to generate a grid of synthetic 
spectra for different metal abundances. To limit the number of 
free parameters, several parameters were frozen. This is the case 
of the stellar mass (15 Mq), the effective temperature at the poles 
*^^eif ~ 34 737 K), the polar radius (6.49 Rq), the projected equa- 
torial rotational velocity (v sin/ - 390 kms '), the inclination 
of the rotation axis (/ - 90°) and the von Zeipel coefficient 
(a - 0.75). The polar temperature and radius are chosen in such 
a way that the average effective temperature and log g of the star 
as seen by the observer coiTespond to 30 918 K and 3.67 respec- 
tively, i.e. agree with the numbers inferred in Sect. 15. 2. Tl This 
grid of models was then used in the same manner as the pure 
plane-parallel models of Sect. 15.2.21 to adjust the three spectral 
regions (see Fig.|6l) by varying the CNO abundances. The result- 
ing CNO abundances are listed in the bottom part of Table|3] As 
can be seen from this table, the differences between the abun- 
dances derived assuming either a spherical star or a rotation- 
ally flattened star are very small and the results of this section 
therefore fully confirm the conclusions of Sect. l5.2.2l This good 
agreement between the CNO abundances stems from the fact 
that in the rotationaiiy flattened model we have chosen the aver- 
age (area weighted) T^ff and log g to match the best-fit parame- 
ters of the pure plane-parallel model. 

Finally, we have tested the sensitivity of the synthetic spectra 
computed from our rotationaiiy flattened models to some of the 
frozen parameters. For this purpose, we have varied the mass to 
8 and 20 Mq, the inclination to 70° and the von Zeipel coefficient 
to 0.50 and 1 .0. In each case, the other parameters were adapted 
in such a way as to recover the mean surface temperature and 
gravity. The resulting differences in the synthetic spectra over 
the three spectral regions adopted for the CNO diagnostics, with 
respect to the best-fit model, were always found to be negligible. 

6. Discussion 

6.1. Couid HD 93521 be a hot subdwarf? 



Ebbets & Savage (T9S2^ suggested that HD 93 521 could be a 
low-mass Population II star. Lennon et al. ( 119911 ) argued that this 
is unlikely because of the strength of the metal lines that would 
be at odds with a low metal abundance expected for a Population 
II star. However, subdwarf O (sdO) and subdwarf B (sdB) stars 
are evolved low-mass objects (post red-giant branch, post hori- 
zontal branch or post asymptotic giant branch) which display a 
wide range of helium abundances and, in some cases, a strong 
enrichment of iron group elements (Heber et al. 120061 Heber 
120091 ). Thus the situation appears less clear cut than previously 
assumed. Therefore, it is worth to critically consider whether 
HD 93521 could be a subdwarf OB star. 

The distinction between sdB and sdO stars stems from the 
presence of He ii lines in the spectra of the latter and their ab- 
sence in the former. Since HD 93521 displays weak but definite 
Hen lines in its spectrum, we have to compare its properties 
mostly with those of sdO stars, although the effective temper- 
ature of HD 93521 that we infer from our model atmosphere fits 
would actually be more typical of an sdB star. 

Let us first consider the abundances. Most sdO stars display 
spectra that point at an He enrichment (up to a factor 10^ with re- 
spect to solar, e.g. Heber et al. 2006 ) which is usually coiTelated 
with enhanced C and/or N abundances (Heber 2009 and refer- 
ences therein). In the spectrum of HD 93521, we find indeed a 
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moderate enhancement of helium and nitrogen (Sect.|5]l, but at a 
less extreme level (for helium) than in typical sdOs. The abun- 
dance pattern of HD 93521 could instead be the result of strong 
rotational mixing in a normal O-type star. Actually, the projected 
rotational velocity of 390 km s"' is much larger than what is ob- 
served in hot subdwarfs. Indeed, with the exception of one ob- 
ject, single sdB stars are slow rotators with v sin / < lOkms"' 
(Geier et al. 2012 ). The exception is EC 22081-1916 which dis- 
plays V sini =163 km s ' and could be the outcome of a merger 
event (see Geier et al. [201 2V Hot subdwarfs in very close bi- 
naries could be spun up by tidal forces, but the binary fraction 
among He-enriched sdOs is rather low (a few percent), much 
lower than in sdBs (Napiwotzki et al. 2004). It seems very un- 
likely also that HD 93521 could be in a close binary system: its 
large v sin / indicates that such a system would be seen under an 
inclination near / ^ 90° and should hence display short-term ra- 
dial velocity variations that would have been clearly detected in 
previous intensive spectroscopic monitoring campaigns (Rauw 
et al. 120081 see also Sect.lOi. 

Hot (Teff > 30000 K) subdwarfs of the sdB category are 
known to display short period (a few minutes) pressure mode 
non-radial pulsations. Much longer periods (45 min to 2 hrs) at- 
tributed to gravity modes were also observed in somewhat cooler 
sdBs (reff < 30000 K; Woudt et al. l2007l see also Fontaine et al. 
12008 and references therein). There is only one known exam- 
ple of a pulsating sdO in the field (Woudt et al. 2007) which 
has very fast pulsations (periods of 1 - 2 minutes). The non- 
radial pulsations detected in HD 93521 have periods of 1.75 and 
2.89 hrs (Rauw et al. 120081 ) much longer than the p-mode oscil- 
lations seen in hot sdBs. Although they are comparable to the 
periods of slowly pulsating sdB stars, the effective temperature 
and gravity place HD 93521 clearly outside the range of this cat- 
egory (see Fig. 1 of Fontaine et al. i2008D . Our model atmosphere 
fits (Sect.|5]) yield an averaged logg - 3.67, much smaller than 
for typical sdOs (logg in the range 4.6 to 6.7, see Fig. 2 of Heber 
et al. 120061 ). 

HD 93521 has an apparent magnitude of V = 7.06. Its 
Hipparcos parallax, (0.85 + 0.49) x 10""^ arcsec (van Leeuwen 
120071 ). indicates a distance that is significantly larger than 
for most other sdOs of comparable apparent magnitude (e.g. 
BD+75° 325, V = 9.44, H = (7.39 + 0.95) x 10"^ arcsec, sd05; 
HD 49798, V = 8.18, H = (1.20 + 0.50) x 10"^ arcsec, sd05.5), 
although the error on the parallax of HD 93521 and HD 49798 
are admittedly quite large. The luminosities of sdOs actually 
span two orders of magnitude (Heber 120091 ) and part of this 
wide range is due to their wide range of temperatures (36000 
- 78000 K, Heber et al. 2006 ). Our determination of the surface 
temperature would place HD 93521 towards the border between 
sdO and sdB stars. This means that we have to compare the par- 
allax of HD 93521 with that of late-type and hence probably in- 
trinsically fainter sdOs, thereby enhancing the disagreement. 

As an additional criterion for the sdO scenario we can con- 
sider the stellar wind features. The winds of sdOs are weak in 
comparison to those of normal O-stars (see e.g. Hamann et al. 
119811 ) and the associated UV spectral features are thus also gen- 
erally weaker. lUE and HST spectra of HD 93521 around the rel- 
evant wind lines are shown by Massa ( 119951 ), Howarth & Reid 
( 119931 ) and Bjorkman et al. (119941 ). There is a rather prominent 
P-Cygni feature in the Civ /1/11548, 1551 doublet which con- 
trasts with the narrow, symmetric absorptions seen in the sdO 
star HD 49798 (Hamann et al.'T98 1). However, the profiles of the 
wind lines (including the Civ doublet) of HD 93521 are rather 
peculiar for O-type stars. These peculiarities have been inter- 
preted as the signature of an axisymmetric wind structure with 



a density contrast of about 60 between the equatorial and polar 
wind (Bjorkman et al. '1994). The determination of the global 
mass-loss rate of HD 93521 is thus not straightforward (spher- 
ically symmetric wind models obviously fail in this context). 
Using their axisymmetric model, Bjorkman et al. (.1994 ) quote 
the best-fit values of the product of the polar mass-loss rate times 
the ionization fraction for C iv. Si iv and N v. These values range 
between 10 and 10 " Mq yr"', which, if the ionization frac- 
tions were close to unity, would be quite comparable to the mass- 
loss rates of the sdOs studied by Hamann et al. ( 198 1). However, 
Bjorkman et al. (119941 ) argue that the ionization fractions are sig- 
nificantly less than unity and that their result would be consistent 
with M - 10"^-^ Mq yr ' . On the other hand, the four sdOs stud- 
ied by Hamann et al. ( 119811 ) are much hotter than HD 93521 and 
a comparison with late sdO or sdB stars might be better. So far, 
observational evidence for winds in this category of subdwarfs is 
scarce and the few cases with a positive detection of wind emis- 
sion in the He line are consistent with mass loss rates of 10 " 
to lO~'''M0yr ' (Unglaub l2008l ). Therefore, our present knowl- 
edge of the wind of HD 93521 does not allow us to draw a firm 
conclusion about the comparison of its mass-loss rate with that 
of hot subdwarfs. 

Finally, we need to address the issue of the X-ray emis- 
sion. Little is known about the X-ray emission of hot subd- 
warfs. An exception are two compact binary systems, HD 49798 
and BD+37° 442, which host a white dwarf or neutron star ac- 
creting material from the sdO component (Tiengo et al. 120041 
Mereghetti et al. | 2009l |20TT , La Palombara et al. 20T2l). The 
X-ray spectra of these objects are rather soft and are well fitted 
by a model consisting of a black-body component (kT x 40 eV) 
along with a power-law (T x 2). Both objects show pulsations of 
the X-ray flux with periods of 13.2 and 19.2 s. As we have shown 
in Sect.m our X-ray data are not well described by this spectral 
model and do not display any signature of pulsations with pe- 
riods between 1.0 and 1000 s. This excludes the possibility that 
HD 93521 could be an sdO with a compact companion. 

During the eclipse of the white dwarf in HD 49798, 
Mereghetti et al. (201 1) detected an X-ray luminosity of 2 x 
lO'^'ergs"', corresponding to Lx/^boi = 7 x 10"** which they 
attribute to the intrinsic X-ray emission of the sdO component. 
Whilst this Lx/Lbo\ ratio is very similar to that of HD 93521, and 
of normal O-type stars in general, overall the various arguments 
presented in this section are definitely at odds with an sdO or 
sdOB nature of HD 93521. 

6.2. The runaway status 

We have measured the radial velocity of HD 93521 on our mean 
ELODIE spectrum. As this spectrum results from the combina- 
tion of many individual spectra taken at various NRP phases, one 
expects it a priori to provide line profiles that should be rather 
free of the NRP signatures. We have only measured lines that 
are expected to be relatively free of strong blends and have a 
rather symmetric shape. Our line list is composed of Hy, Hj3, 
He I AA4471, 4922 and Hen AA420Q, 4686, 5412. In this way, 
we infer a heliocentric radial velocity of 9.5 + 6.3kms"', which, 
using the formaUsm of Moffat et al. (1998) and Naze (2004]), 
translates into a peculiar radial velocity of 10.0 ±6.3 km s"', well 
below the conventional 30kms ' threshold (Cruz-Gonzalez et 
al. 1 19741 ) for runaway stars. 

Using radial velocities compiled from the literature, Gies 
(1987) found a heliocentric radial velocity of -13.9kms"' cor- 
responding to a peculiar radial velocity of -12.5 km s"'. Whilst 
this result is still below the threshold for a runaway status, it 
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is actually quite different from our result. This prompted us 
to check the RVs on the average Hei line profiles in the data 
used for the analysis of the NRPs (Rauw et al. 2008 ). The lat- 
ter yield heliocentric RVs of -16.2 and -2.9kms ' for the Hei 
/1 5876 line respectively in February 2006 and April 2007, as 
well as 35.6 and 3.8 km s ' for the Hei /1 6678 line respectively 
in February 1997 and May 2005. Finally, measuring the same 
lines as for the mean ELODIE spectrum on our SOPHIE spec- 
trum yields a heliocentric RV of 21. 4+11.1 km s" ' . Whilst these 
differences could be due to orbital motion in a (most likely) long- 
period (years) binary system, an alternative explanation could 
be the NRPs. Rzaev & Panchuk (2006) studied the bisector 
radial velocities of a number of Hei lines in the spectrum of 
HD 93521 over a time span of 2.7 hrs. They found apparent vari- 
ations of the radial velocities with amplitudes of tens of kms ^ 
Therefore, the radial velocities of HD 93521 could strongly de- 
pend on the specific phase of the NRPs or the sampling of the 
latter. This means that the true heliocentric RV is difficult to as- 
sess. Nonetheless, in view of the above results, it seems unlikely 
that the true peculiar velocity of HD 93521 could be significantly 
larger than the 30 km s ' threshold. 

The Hipparcos catalogue based on the revised reduction (van 
Leeuwenl^UOT^ lists a pai-allax of H = (0.85 + 049) x 10"^^ arcsec 
and proper motions of /i„ cos (5 - (0.32 + 0.40) x lO^^* and 
lis - (2.44 + 0.37) X lO^^* arcsec yr"'. We have used these num- 
bers to estimate a peculiar tangential velocity of HD 93521 of 
V, = 26.8 ± 9.2 km s"'. Adopting the criterion of Moffat et al. 
(I1998I I. this number is well below the threshold for considering 
HD 93521 as a runaway star. 

Based on its kinematic properties, HD 93521 can thus not 
unambiguously be classified as a runaway star. However, de Wit 
et al. ( 120051) argue that stars located more than 500 pc above 
the Galactic plane must de facto be runaways, as there is cur- 
rently no evidence for star formation in the halo. Indeed, in 
their study of a sample of 10 high Galactic latitude B-type stars 
with masses between 5 and 13 M©, Ramspeck et al. (j2001j did 
not find a conclusive candidate for a star formed in the halo. 
Although they were lacking kinematic data for some of their 
targets, they found that the evolutionary time scales were suf- 
ficiently long for the stars to have travelled to their current loca- 
tion if their ejection velocity reaches values up to 440 km s"' . To 
estimate the evolutionary age of HD 93521, we have plotted the 
star in a (log T^^, \ogg) diagram against the evolutionary tracks 
and isochrones of Ekstrom et al. (2012). The result is shown in 
Fig.[8] HD 93521 falls between the evolutionary tracks of initial 
mass 20 and 25 Mq and close to the 6 Myr isochrone. Note that 
this estimate might actually be a lower limit of the true age, if the 
evolution of the star was affected by binary effects, such as mass 
transfer during a Roche lobe overflow phase or a merger event 
(follo wing ejection through dynamical interactions in a dense 
cluster 'FI leading to a blue straggler star (de Mink et al.'2010'). 

Dynamical interactions in dense clusters can produce very 
high maximum ejection velocities of up to 700 km s"' for a 
60 Mq O-star and even up to 1400kms"' for low-mass stars 
(Leonard I1991I I. whilst the supernova ejection scenario gen- 
erally predicts maximum ejection velocities of < 300 km s"' 
(Portegies-Zwart 2000). However, Przybilla et al. (2008) pro- 
posed that an ejection velocity of about 400km s"' could re- 
sult from a supernova explosion in a tight Wolf-Rayet + B-star 
binary. From a sample of 96 likely OB runaway stars, Silva 



In the dynamical interactions scenario, the ejection of a close bi- 
nary system is a rare event and should proceed at a comparatively low 
ejection velocity (Leonard & Duncan ri990b . 



3 



3.5 



4 









/ / ^/ / / / / ^ 


/ 


/ // / / / / A 


/ ; / 


/yA'^r' / / / / 

/ /i/y'/jP /^/ // - 


~ / ', / . 

/ '■''A 

/ ^ h '' 

/ / / y 




/ / / ; / 

J / j ^ p 

h r A 
// r / 

_ // n 1 


i ' /' " / '' ''/ ' 


32 ' 25 ' 20 'U /P // 

12 9^7 



4.6 4.5 4.4 4.3 

l°e T^„ (K) 



Fig. 8. Location of HD 93521 in the (log Teif, logg) diagram 
(blue dot with error bars) as inferred from our model atmosphere 
fit (Sect. l5.2Tn The main-sequence parts of the Geneva evolu- 
tionary tracks from Ekstrom et al. ( I2012I I are shown by the solid 
lines for stars initially rotating at 40% of the critical velocity. 
The initial mass of the star is given in Mq by the label at the bot- 
tom of the track. The dashed lines yield the same tracks for zero 
rotational velocity. The isochrones for ages of 6.3 and 7.9 Myr 
are shown by the dotted red fines. 

& Napiwotzki (1201 11 1 found evidence for a bi-modal distribu- 
tion of tfie ejection velocities witfi two populations witfi veloc- 
ities < 300kms"' and around 400 - 500kms"' respectively. 
Tfierefore, assuming an ejection velocity of 400 km s"', tfie min- 
imum time requested to reacfi tfie star's present location, 1 kpc 
above tfie Galactic plane, would be about 2.4 Myr, wfiicfi is well 
below tfie estimated evolutionary age. Tfierefore, in ligfit of our 
current results, we cannot rule out tfie possibility tfiat HD 93521 
could be a runaway star, probably observed near tfie apex of its 
orbit. 

7. Summary and conclusions 

Our analysis of tfie X-ray and optical data of HD 93521 provides 
new clues about tfie nature and origin of tfiis star. We found no 
evidence of a stellar cluster or otfier tracers of star formation 
activity wfiicfi makes it unlikely tfiat tfie star fias formed at its 
current position. Tfie X-ray data did not reveal any evidence for 
a compact companion tfiat could fiave been produced in a su- 
pernova explosion. Tfie stellar parameters inferred from an anal- 
ysis of tfie optical and X-ray spectrum are consistent witfi an 
'ordinary' Population I late O-type star and reject an interpre- 
tation of tfie star as a low-mass subdwarf. Tfie He and N sur- 
face abundances are enfianced, as expected for a fast rotating 
evolved main-sequence O-type star (Meynet & Maeder 200()]l. 
Tfie X-ray spectra can be represented by a two-temperature op- 
tically tfiin tfiermal plasma model witfi a dominant component 
witfi kT a; 0.3 keV and a ratfier fiot secondary plasma component 
witfi kT a; 3 keV. At first sigfit, tfie subsolar metallicity inferred 
from fitting tfie X-ray spectrum could be somewfiat of a surprise. 
However, it is not uncommon to derive subsolar metallicities 
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from the fit of CCD X-ray spectra of O-type stars. A prominent 
example of this situation is the population of OB stars in M 17 
(Broos et al. 2007 ) which has a metallicity of x 0.3 Zq. 

Whilst the runaway status of the star remains ambiguous, 
the most likely explanation seems nevertheless that the star has 
not formed in the Galactic halo, but was rather ejected from the 
plane, either through dynamical interactions in a dense cluster 
or an unstable multiple system, or as a result of a supernova ex- 
plosion in a binary system. The rapid rotation of the star, if not 
due to mass and angular momentum transfer in a binary system, 
could be the result of the dynamical interaction, maybe associ- 
ated with a merger event (e.g. de Mink et al. 2010) that would 
produce a blue straggler and could alleviate the flight time con- 
straint. One open issue that remains with this scenario is how 
the star maintained its extreme rotational velocity well after the 
event that led to its ejection from the plane. A better knowledge 
of the proper motion and distance of the star, as will be provided 
by ESA's forthcoming GAIA mission will certainly be extremely 
beneficial to complete our understanding of this intriguing star. 
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